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The ac Josephson effect in a ferromagnetic Josephson junction, which is composed of two su- 
perconductors separated by a ferromagnetic metal (FM), is studied by a tunneling Hamiltonian 
and Green’s function method. We obtain two types of superconducting phase dependent cur- 
rent, i.e., Josephson current and quasiparticle-pair-interference current (QPIC). These currents 
change their signs with thickness of the FM layer due to the 0-7T transition characteristic to 
the ferromagnetic Josephson junction. As a function of applied voltage, the Josephson critical 
current shows a logarithmic divergence called the Riedel peak at the gap voltage, while the 
QPIC shows a discontinuous jump. The Riedel peak reverses due to the 0-7T transition and 
disappears near the 0-7r transition point. The discontinuous jump in the QPIC also represents 
similar behaviors to the Riedel peak. These results are in contrast to the conventional ones. 

KEYWORDS: ac Josephson effect, Riedel peak, ferromagnetic Josephson junction, superconductor, ferro- 
magnet, 0-state, pi-state 



1. Introduction 

The Josephson effect is a macroscopic quantum phe- 
nomenon involving phase coherence between two super- 
conductors (SC’s), and is characterized by a zero voltage 
current through a thin insulating barrier (dc Josephson 
effect). On the other hand, when a finite voltage is ap- 
plied across the junction, an alternating current flows ac- 
cording to time-dependence of phase-difference, 6, in two 
superconductors (ac Josephson effect), which gives two 
types of current as, /ji = / c isin0 and Ij 2 = I^cosO. 1 
The former is the same as the dc Josephson current ex- 
cept for time-dependence of 0, while the latter is a phase- 
dependent dissipative current inherent to the ac case. 
As a function of applied voltage, V, the amplitudes, I c i 
and I c 2 , show singularities originating from the gap, A, 
in an s-wave superconductor. The I c \ has a logarithmic 
divergence called the Riedel peak at the gap voltage, 
V = 2 A/e, due to the square root singularities in the 
density of states , 2-7 while the I c 2 shows a discontinuity 
at V = 2 A/e and is zero below V < 2 A/e at zero tem- 
perature . 3,5 The I c 2 is experimentally observed near the 
superconducting transition temperature, T c , 8 in various 
types of Josephson junctions . 9-1 4 

Recently, a Josephson junction with a ferromagnetic 
metal (FM), which is called the ferromagnetic Josephson 
junction, has been actively studied both experimentally 
and theoretically . 15-28 One of the most interesting phe- 
nomena in the ferromagnetic Josephson junction is the 
oscillation of the I c 1 as a function of the thickness of 
the ferromagnetic film. The mechanism of the oscillation 
is similar to that of Fulude-Ferrell-Larkin-Ovchinnikov 
(FFLO) state . 29,30 Cooper pairs penetrating into the FM 
acquire a finite center of mass momentum proportional 
to the magnetic exchange splitting, h ex , between up- and 
down-spin bands. As a result, the pair correlation in FM 
oscillates as a function of the thickness of FM. If the 
thickness of the FM is about a half of the period of 
the oscillation, the current-phase relation is shifted by 



7 T from that of a conventional Josephson junction (0- 
junction) like as, I c 1 < 0. This is called a 7 r-junction, 
which has potential applications as a quantum bit . 31,32 
In addition, the 7 r-junction is used in some experiments 
to measure a nonsinusoidal current-phase relation. Using 
the ac Josephson effect, Sellier et al. experimentally evi- 
denced the second harmonic term given by sin 26 in the 
SC/FM/SC junction . 20 However, most of studies on the 
7 r junctions have been so far focused on the dc Joseph- 
son effect. Studies on the ac Josephson effect in ferro- 
magnetic 7 r junctions will open a new pathway of basic 
physics and will contribute to realize the quantum bit 
including SC/FM/SC junctions. 

In this paper, we study the ac Josephson effect in a 
SC/X/SC junction, where X is either a NM or a FM. 
Using a tunneling Hamiltonian and Green’s function 
method , 33-35 we obtain two types of phase dependent 
current, i.e. the ac Josephson current (Iji = I c isin$) 
and the QPIC (Ij 2 = / C 2 cos 0 ). In a SC/NM/SC junc- 
tion, 7 c i and I c 2 monotonously decreases with the thick- 
ness of NM. As a function of applied voltage, I c 1 shows 
a logarithmic divergence (the Riedel peak) at the gap 
voltage, V = 2 A/e, while I c 2 discontinuously jumps at 
the same voltage. On the other hand, in a SC/FM/SC 
junction, I c 1 and I c 2 exhibit the strong dependence on 
the thickness of FM and changes their signs by crossing 
the 0 - 7 T transition point due to the magnetic exchange 
splitting between the up- and down-spin bands in FM. 
In particular, it is predicted that the Riedel peak in I c 1 
disappears at the 0 - 7 T transition and the 0 - 7 T transition 
occurs in I c 2 like as I c ±. 

The rest of this paper is organized as follows. In §2, 
we introduce the model Hamiltonian including the tun- 
neling Hamiltonian, and explain the formulation to cal- 
culate the ac current by the thermal Green’s function 
method. In §3, the ac Josephson current and the QPIC 
in the SC/NM/SC and SC/FM/SC junctions are shown 
as functions of thickness and applied voltage. We com- 
pare these two types of junctions and discuss similarities 
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and differences. Summary and discussion are given in §4. 
Below, h = 1 and fce = 1 are used in the equations. 



2. Tunneling Hamiltonian approach for ac 
Josephson current 



2.1 Model Hamiltonian 

We consider a junction composed of a ballistic FM 
with the thickness d and s-wave SC electrodes as shown 
in Fig. 1. A dc bias voltage, V, is applied across the 
junction. 

The SCs and the FM are connected by tunneling 
Hamiltonian. The total Hamiltonian is given by 
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Fig. 1. Fourth order diagram with a tunneling matrix element 
contributing to Ij. Insulating barrier is at interfaces between SC 
and X. 



where (J) is the expectation value of the current opera- 
tor. The first term of eq. (7) is the Josephson current 
and / c i is the Josephson critical current. The second 
term is called the quasiparticle-pair-interference current 
(QPIC). 

2.2 Josephson critical current and QPIC formula 
In the SC/FM/SC junction, the fourth order term of 
J with regard to T r y is shown in Fig. 1. Detailed cal- 
culations are given in Appendices A and B. For F ex //i, 
bJn/ H « 1, and temperature, T=0 [K], I c \ and I c 2 are 
given by, 

G(Ao-\E-eV\) 

^E 2 -A^A 2 ~(E-eV) 2 



Ir 1 = 



e>o A, 



dE 



' Aq 



+ h.c., (5) 

where ipi tl7 = (r) is the electron field operator with 

the position r in the region i (=L, R, or FM) and the 
spin a (=|,J,). We adopt the BCS mean field Hamilto- 
nian with the s-wave gap A and the phase variable 

$L(R) in the left (right) SC. Then, the phase difference is 
given by 6 = (?l — 0r. The electron mass and the chemi- 
cal potential are denoted by m and /i, respectively. Note 
that =FeT/2 is added in the Hamiltonian (2), since the 
applied voltage, V , imposes a chemical potential differ- 
ence, V/2, at each boundary between FM and SC. The 
Hamiltonian of the FM, Hfm, which has no impurity 
scattering, has the exchange energy, h ex . Ht is the tun- 
neling Hamiltonian, whose matrix element is denoted by 
T r/r i and has a finite value at the SC/FM boundary as, 
T r y — To<5 (r — r') 6 (r — r L ( R )). t'l(r) is the position 
of the interface between the left (right) SC and the FM. 

We calculate the expectation value of a current oper- 
ator, 



J = — ie^ J drdr'T r ye leVt ^ 2 ipF^'ifFM,cr + h.c, 

(6) 
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where x involves both time, t, and r. The superconduct- 
ing phase dependent current is given by, 

I j - (^) = hi sind + I c2 cos 6, (7) 
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where Ao is the superconducting gap at T=0 [K], <jq = 
167re 2 ToIVL(0)IVR(0) [mv/{2nd)\ i is a constant deter- 
mined by materials and interface, AA(r) is the density of 
states in the left (right) lead at the Fermi level, up is the 
Fermi velocity, v and d are the volume and the thickness 
of the X, respectively. Ci(x) and Si(a:) are the cosine and 
sine integrals, respectively. When h ex = 0, our formula- 
tion reproduces the current in the SC/NM/SC junction. 

For d/f 0 >> 1, eq. (8) is simplified as 
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Fig. 2. (a) Normalized Josephson critical current, 7 c i, and (b) the 

quasiparticle-pair-interference current, I c 2 , as functions of d/£o 
for a SC/NM/SC junction ( h ex = 0), where / c i = I c ine/croAo 
and I C 2 = / c 27re/aoAo. 




Fig. 3. (a) Normalized Josephson critical current I c i and 

(b) quasiparticle-pair-interference current I c 2 as functions of 
eV / 2 Ao in a SC/NM/SC junction , where I c 1 = fciW^oAo 
and I C 2 = / c 27re/croAo. 
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where £o = Uf/ 2A 0 . It is found that I c i decays with d in 
the power law as 1 / d. This behavior is consistent with the 
case of dc Josephson critical current in a double barrier 
Josephson junction. 28 On the other hand, for d — > 0, eqs. 
(8) and (9) reproduce the current in SC/I/SC junctions. 5 
The integration in eqs. (8) and (9) is numerically carried 
out and results are shown in the next section. 



3. d- and R-dependence of J c i and J c 2 

In this chapter, we show thickness d and bias voltage V 
dependences of 7 c i and I c 2 . First, the SC/NM/SC junc- 
tion is discussed to see the difference between NM and 
FM cases clearly. Then, we will give detailed discussions 
on the SC/FM/SC junction. 



3.1 SC/NM/SC junction 

The thickness dependence of I c 1 and I c 2 is shown for 
several values of V in Fig. 2. In Fig. 2 (a), the vertical 
axis is the normalized I c 1 and the horizontal axis is the 
thickness d normalized by the coherence length £q- It is 
found that I c 1 shows a monotonic decrease as a function 
of d. For d/£o » 1, I c i decreases in the power law as 1 /d 
shown in eq. (10). In Fig. 2 (b), I c 2 is plotted. It is found 
that I c 2 decreases as a function of d for eV/2A a = 1.5, 
while it is always equal to zero at eV / 2 Aq = 0. 



The 7 C 1 in the SC/NM/SC junction is shown for sev- 
eral values of d in Fig. 3 (a). The vertical axis is the 
normalized I c i and the horizontal axis is the normalized 
voltage, eV/2A 0 . In Fig. 3 (a), 7 c i shows the Riedel peak 
at the gap voltage similar to that in a SC/I/SC junction. 
In this system, the Riedel peak exhibits weak dependence 
on d. The I c 2 is shown in Fig. 3 (b). It is found that J c 2 
has discontinuity at eV / 2 A$ = 1. The behavior of I c 2 is 
the same as the case of SC/I/SC and 7 C 2 exhibits very 
weak dependence on d. 

Here, we examine the above results in details. 7 c i shows 
the monotonic decrease as a function of d , which repre- 
sents a decoherence of the Cooper pair penetrating into 
the NM. This behavior is equivalent to the case of dc 
Josephson effect. On the other hand, the behavior of 
7 C 2 is quite different from that of I c i, since I C 2 has a 
finite value only when V is larger than the gap volt- 
age as shown by Fig. 3 (b). Therefore, it is interpreted 
as quasiparticles in the band below the superconduct- 
ing gap contribute to I C 2 and carry the superconduct- 
ing phase coherent dissipative current as in the case of a 
SC/I/SC junction. 36 Moreover, I C 2 in a long SC/NM/SC 
junction (d » £ 0 ) oscillates as a function of d and V 
with the period ( AnE/v^d because of cos in eq. 

(9). This behavior represents the phase shift by 7r and 
the 0-7T transition occurs as a function of d and V. In 
the R-dependence of I cl , the divergence of 7 cl at the gap 
voltage is also shown in the SC/I/SC junctions. Since, the 
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Fig. 4. (a) Normalized Josephson critical current, I c l, and (b) 

quasiparticle-pair-interference current, I c 2 , as functions of d/£o 
for h ex /A o — 50 in SC/FM/SC junction, where J c i = 
/ c i7re/croAo and I c 2 = 7 c 27re/croAo. 
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Fig. 5. (a) Normalized Josephson critical current, I c 1 , and 

(b) quasiparticle-pair-interference current, I c 2 , as functions of 
eV/2Ao for h ex /Ao = 50 in SC/FM/SC junction, where f r \ = 
/ c l7re/(ToAo and / c2 = I C 2 ne/a 0 A 0 . 



density of states has a square-root singularity at the gap 
edge as we can see in eq. (8), the amplitude of the Joseph- 
son critical current produces a logarithmic singularity at 
the gap voltage. This mechanism on the singularity is 
the same as that of SC/I/SC junction. In J c2 , the current 
vanishes for eV/2A(, < 1 and jumps at the eV/2Ag = 1. 
This behavior is similar to that of the SC/I/SC junctions. 
In a long SC/NM/SC junction (d » £ 0 ), / C 2 oscillates 
as a function of V with the period (^itE /vp)d as shown 
in eq. (9). 

3.2 SC/FM/SC junction 

In this subsection, the case of the SC/FM/SC junc- 
tion is compared with that of the SC/NM/SC junction 
from the viewpoints of the (/-dependences of I c 1 and I c 2 
as shown in Fig. 4. In Fig. 4 (a), the normalized I c 1 is 
plotted as a function of d/£ o- For eV/2A 0 = 0, i.e., dc 
Josephson effect, I c \ shows a damped oscillatory behav- 
ior as a function of d and the 0-7T transition occurs by 
increasing d unlike the case of SC/NM/SC junction. For 
eV/2A 0 = 1.5, / c i is qualitatively the same as that in the 
dc Josephson effect. On the other hand, / c2 is zero for 
eV/2A 0 = 0, but has a finite value for eV/2A 0 = 1.5 
similar to the SC/NM/SC junction. However, the / c2 
changes its sign with increasing d due to the 0-7T tran- 
sition and vanishes at the transition point as shown in 
Fig. 4 (b). 

Figure 5 shows the ac current amplitude as a func- 
tion of V in the SC/FM/SC junction. In Fig. 5 (a), the 



vertical axis is the normalized I c 1 and the horizontal 
axis is the normalized voltage, eV/2Ag. In this system, 
the Riedel peak exhibits a strong dependence on d and 
changes its sign with increasing d due to the 0-7T tran- 
sition. Therefore, near the thickness at which the 0-7T 
transition occurs, the Riedel peak disappears as shown 
in Fig. 5(a). In Fig. 5 (b), it is found that / c2 has a finite 
value above eV/2A 0 = 1. For d/£ 0 = 0.01, the behavior 
of / c2 is the same as those of SC/I/SC and SC/NM/SC 
junctions. On the other hand, / c2 changes its sign with 
increasing d due to the 0-7T transition and disappears 
near the thickness at which the 0-7T transition occurs, as 
shown in Fig. 5 (b). 

The I c i shows the damped oscillatory behavior as a 
function of d and the 0-7T transition occurs. In eq. (8), 
only the ratio of h ex and vp determines the period of os- 
cillation in the / C i-(Z curve. The mechanism of 0-7T tran- 
sition in the ac Josephson effect is the same as that of 
the 0-7T transition in the dc Josephson effect. Concerning 
/ c2 , it has a finite value above eV/2Ao = 1, similar to 
the SC/I/SC and SC/NM/SC junctions. The behavior of 
I c 2 is quite different from that of SC/NM/SC junction. In 
the SC/FM/SC junction, the oscillating term in eq. (9) is 
composed of two part. One is cos ( 2Ed/vp ) (the region of 
E being from 0 to eV / 2— A 0 ), the other is cos ( 2h ex d/vp ). 
In the practical case, h ex is experimentally larger than 
eV . 17-27 Therefore, the period is dominated by h ex and 
is short compared with £o- The H-dependence of I c \ and 
/ c2 in the SC/FM/SC junction shows remarkable plre- 
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Fig. 6. Ac current I as a function of t/ Ao for /i ex /Ao = 100, 
d/S o = 0.07, eV/Ao = 1.5 and 2.5. / = Iire/ao/S.^ 



nomenon. In I c i, the Riedel peak exhibits a strong de- 
pendence on d. It changes its sign with changing d and 
disappears at the 0-7T transition point as shown in Fig. 
6(a). I C 2 also changes its sign with changing d and disap- 
pears at the 0-7T transition point as shown in Fig. 6(b). 
We expect that these results provides a new method to 
observe the 0-7T transition in SC/FM/SC junctions. 

We take a look at the total ac current given by, 

1 = \l I 2 C 1 + J c 2 sin (2 eVt + x) (11) 

X = arctan (/ c2 // c l) (12) 

Fig. (6) shows the current-phase relation(CPR) for 
h ex / A 0 = 100, d/^ = 0.07, and eV/A 0 = 1.5 and 2.5. 
For eV/2A 0 < 1, only the Josephson current flows be- 
cause of J c2 = 0. Therefore, the CPR represents a con- 
ventional behavior as in a SC/I/SC junction as shown 
by the dashed line in Fig. (6). On the other hand, for 
eV/2A 0 > l,the phase of CPR is shifted by x due to 
the finite QPIC. The behavior is shown in the solid line 
in Fig. (6). From these results, we can confirm the exis- 
tence of the QPIC when the voltage dependence of the 
current-phase relation is measured. 

4. Summary 

In summary, we have studied the ac Josephson ef- 
fect in the SC/X/SC junction (X = NM, FM). Using 
a tunneling Hamiltonian and Green’s function method, 
we obtained the Josephson current and the quasiparticle- 
pair-interference current (QPIC). The Josephson critical 
current, I c i, shows the Riedel peak at the gap voltage, 
V = 2A/e. In the SC/NM/SC junction, the Riedel peak 
exhibits weak dependence on d. The amplitude of QPIC, 
/ C 2 , has a finite value above the gap voltage at T = 0[K], 
These behaviors are similar to those of SC/I/SC junc- 
tions. On the other hand, the critical currents in the 
SC/FM/SC junction show quite different behaviors com- 
pared with those in the SC/I/SC and SC/NM/SC junc- 
tions. In I c i, the Riedel peak exhibits a strong depen- 
dence on d and changes its sign with increasing d due to 
the 0-7T transition. Therefore, at the thickness of the 0-7T 
transition, the Riedel peak disappears. I C 2 also shows the 



0-7T transition with increasing d and vanishes at the thick- 
ness of the 0-7T transition. The ac Josephson current has 
a strong dependence on the applied voltage, and at the 
gap voltage the amplitude of / c i shows logarithmic diver- 
gence. The amplitude of the higher harmonic Josephson 
current also show the applied voltage dependence, similar 
to the Josephson current. The study of the ac Josephson 
effect in the SC/FM/SC junction gives a possibility to 
observe the higher harmonic Josephson current. 

We have studied the ferromagnetic Josephson junction 
in the clean limit in this paper. In many experimental 
situations, FMs are usually in a diffusive transport re- 
gion. However, the essence of 0-7T transition in the ac 
Josephson current may be the same even if the case of 
a dirty FM. Recently, the current-voltage characteristic 
was studied by using the superconducting phase differ- 
ence coupled with the spin wave excitation in the ferro- 
magnetic Josephson junction within the phenomenologi- 
cal model. 37 And the dc Josephson current coupled with 
the spin wave excitation was also discussed. 38 Therefore, 
it might be important to consider the spin wave excita- 
tion for the ac Josephson current from microscopic view 
point. These problems will be left in a future issue. 
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Appendix A: Perturbative calculation of I c i and 

I c 2 in path integral formulation 

A.l Basic formula 

In the path integral formulation, the partition function 
is given by 

Z = Tr[e _/3jF/ ] = / V^*V^e~ s , 



J dr (^J drip* d T tp + H (r)^ , 

(A-l) 

where Wick’s rotation is performed, it — ► r, eU/2 — > 

(/ 3 = T _1 , r is an imaginary time variable, H is a tempo- 
rary boson’s Matsubara frequency. That becomes eV/2 
by an analytic continuation.) and = iffd-rifi, + 

4’FM^r'f’FM + It is noted that the energy of the 

left-hand superconductor (L) is different from that of the 
right-hand superconductor (R), since the voltage, V, is 
applied. When these systems are separated each other, 
the time-development in the left-hand side is given by, 

Ul C t , t') = exp [-* (H L - A L ) ( t - t') 

+ i (eV/2) N l (i - If)] , (A-2) 

while that in the right-hand side is 

Ur (t, t') = exp [-* (Hr - Ar) ( t - t') 

-i(eV/2)N R (t-t')\, (A-3) 



= / P'P*X>'Fexp 
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Here, we make a transformation as, 

tph (r, t) = e*( eV / 2 ) t ^ L (r, t ) , 

V’R I 71 , t) = e _i(eV/2) Vr (r, f) . 
Therefore, H in the action, S, is transformed as 

^bcs = J ~ ^ V 2 A.cr 

-PL J dx^L,T^L,iV’L,J,^L,T> 

ffFM = (unchanged) , 
ffncs = (L — > R) , 



The inverse of the Green’s function for electron system, 
G _ , is given by, 

'-GZ f 0 

G — rp* A— rp 

1 ^FM 1 ’ 

0 T* -Gn 



— “^0 + 
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+ J dxdx'T X:X fe ieVt/ 2 ^y FM ,a 

a reR.r'eFM 

+ h.c. (A-9) 

For convenience, the tilde is abbreviate. Where ip\ J. = 

(®), i is L, R, or FM, x = (r,r), and T x y = 
T r y5(r — t'). The current operator is defined as, 

J = — = — e [dr^-ie\[drn F ,H , 
dt I dt I 



' 0 T 0 
P= T* 0 f 
_ 0 T* 0 

-Gr = kd 0 - ( 



(A-19) 

-A + <T 3 - A L e if?L<T3 <7i 
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GpM d T d o A t t iah ex <73 



=~ ie Yl 
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+/i.c., 
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(A-20) 

(A-21) 



xS (x — x') , (A-22) 

» " p hifir-p O 

T= n _ pi Q. T * 5(r-r'). (A-23) 

g 1 , 173 are the Pauli matrices and do is a unit matrix. 
6 * L (r, is the phase of superconducting order parameter, 
A L (R)e ieL(R) . Below, we consider the case that the super- 
conductor L and R are same. 

An auxiliary held, <j>, that couples to the current oper- 



(A-10) 

An auxiliary held, q b, that couples to 
(A- 11) a t° r i J> is introduced in ecp (A- 10) as, 



We eliminate the quartic interaction term in Hamil- ‘~’' cond J dx ^ + j ax ^ 

tonian eqs. (A-6) and (A-8) using Stratonovich-Hubbard 

transformation, + f dr 0 J} (A-24) 

r fi/* 1 Jo 

1= DA* DA exp — / dx (A — giptipt) (A* — giptipO , r r , 

J L dJ n ™ T i; J <f>J = (j) dxdx'ty* (x) J H/ (x') , (A-25) 

(A-12) ^ 

where (f> = (r) and J = J (x, a:'). The current operator is 
where A ( *i = A^*\x). With this transformation, the ac- denoted as 
tion becomes, 

0 —ieR 0 

S = S cond + Seh (A-13) J= i e R* 0 0 , (A-26) 
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' tr ' r ' in “ S(t-t'). (A-27) 

S e i = / d® / d® 7 ** (®) [G“] (® 7 ) . (A-15) L r ’ 7 '' J 

^ Tracing out the electron helds from ecp (A-25), we obtain 

The Green’s function, G~ = G~ (x, x 1 ), is a 6 x 6 matrix the effective free energy, F e s, with Josephson junction as 

spanned both in the space of L, FM, R and Nambu space. following, 

The electronic helds are expressed as, r r - a 3 

, . 0F eff = -Tr \ In G" 1 + <j>J 

ftp LT \ t L J J 

= — Tr { In [ (-Go + (p J] + p] ) - (A-28) 

$(3)= ( A . 16) L LV 1 JJ 

Z’fmj. ’ v 7 In eq. (A-28), Tr means taking trace with respect to the 

"Zri matrix element of the Green’s function and integrating 

VZ’ri / with respect to x, x' . In the fourth order perturbation 



(A-16) 



-Tr{ln[(-Go + 0j) + p] } . (A-28) 
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theory about the tunneling matrix, the currei 
proximated to be, 



(J) = 



dF e g 



dcj) 



~ — Tr 

P 



«.= o 



JGri 






+ /? Tr 



JGnPGr 



+^ Tr 



JG 0 PG 0 PG 0 PGo 



The first and second terms doesn’t contribute the 
son current, whose leading term is the third oi 
(A-29 ). 

Because Tr is not affected by a change of basis, 
trix element can be transformed from (r,r) comp 
its Fourier component Here, iu> n is the 

ara frequency of the electrons. Here we made u; 
fact that the tunneling matrix element differs fro 
for r ~ r'ancl in the neighborhood of the SC/FA 
ary, i.e., T r y = T 0 S (r - r') S ( r - r L(R) ), whe 
is the position of the interface between the lef 
SC and the FM. Tracing out the matrix elemei 
Green’s function, we can obtain the current foi 
following, 



h 





^FM ) r 



X /l (&L? i^n + At) ^ ^^FM’ i^n'j 

X 0FM,X ^FM) /r (&Ri i^n ~ in) ^A-30) 

where r = r R — r l and each Green’s function is given by 

A L e i¥,L 



fh(kL,ico n +ifl) = 



(ui n + H) + ££ + A 2 



(A-31) 



/ R (fc L , ico n - in) = (L - R) , (A-32) 

ffFM,t (^FMf ) * w n) = 7 > (A-33) 

l^n ~ SFM 

5fm,x( — ^fm|) — iu n ) = 7 — r , (A-34) 

VjJn + SFM 

a ( R ) = ^-M, (A-35) 

?FM = - A* - CT/lex- (A-36) 



Note that /cfmt and &fmx are independent of each other. 
After integrating Green’s functions of the FM as to 
1 :fm,it, we can obtain the following formula, 



/j = 16eT c 



mv 



1 2 



0 L27TdJ 
1 



Im 



.A 2 

’ P 



E 

fcL,fcR,o;n>0 

1 



K + ^) 2 + E 2 kL K - n ) 2 + E 2 kR 
sin 6 



x cos 



( ‘Zhex \ e ~2cj n d/v F 

V ^'F d J 




Fig. A-l. Paths in the complex plane for the contour integration 
of eq. (Bl) 



16eTn 



mv 



1 2 



L 27tcZ J 
1 



Im 



A E 

^L 5 ^R^n>0 

1 



(Ldn 



X COS 



( 



-nf 

2 h ex 



El(u; n -nr + El R 



-2tJ n d/vF 



\ vp d 
= / c i sin 6 + J c2 cos 6, 



cost 



(A-37) 



where v is a volume of the FM, k F is the Fermi wave 
number, (3 = T _1 . 

Appendix B: Summation of Matsubara fre- 

quency 

We evaluate Matsubara frequency summation in 
(A-37). This summation is performed by a contour in- 
tegration in Fig. B. The integral has the form 

dz 

ic 2 ni J 

1 



Ic = ~ 

m = 



£^-f(z)M z ) ei2zd/v *, (B-l) 



1 



z T %n — E k L z T in T E k L 

1 1 


z — in — E kR z — in + E kR 
C = C 1 + C 2 + C U , + C R , 



(B-2) 
(B-3) 

where n F (z) is a Fermi distribution function, Ci and C 2 
are a closed path enclosing the pole of the anomalous 
Green’s function, C w is a integration on the real axis of 
complex plane, and C R is a large semicircle of radius R 
in the limit as R — > oo. The contour integration with 
Cr is zero because of the Jordan’s theorem. When each 
contour integration is performed , the following result is 
obtained, 

Ici= 1 



^E k L E kn 



1 



1 



i2n + E kR — E kR 

i2n 



x n F (E kR ) exp ( i 
1 



i2n + 

2 E kR 



Ek K ~ Ej 



k-L 



v F 



(B-4) 



/C2 = - 



4 E klj E kR 
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Author Name 



1 



i2Q — E kR — E kh i2Q — E kR + E kh 
(,i2Q-2E kR \ 

x n F {-E kK ) exp i d ) , 



fF 



(B-5) 



L, = 



dm 



1 



1 



I —oo 22^2 -£//cr 

l l 



i2Q + 2m - 2E kh i2Q - 2m - 2 E kR 

■°° dm 1 1 

-oo ^^k R E kR z 2(1 } E kR E k L 

1 1 



z2fi + 2w - 2S fcL z2fl - 2w + 2E kR 

■°° 1 

-oo ^ ^■^ki J E kR z 2(1 E kR T E k L 

1 1 



z2(2 -|- 2uz -t~ 2 E kR z2fl — 2m — 2 E kR 

■°° 1 

-oo ^^k R E kR z 2(1 -(- E kR T E k L 

1 1 



z2(2 + 2w + 2S fcL z2fl - 2m + 2 E, 



k R J J 



xrip (w) exp ( i—d I . 

V F 



(B-6) 



Here, we consider the case of absolute zero temperature. 
Then zip (.Efc L(R) ) = 0, n F (m) = 0(— m). Making an ana- 
lytic continuation z2(7 — > eV + iS and integrations with 
respect to m and E kh or E kR1 we obtain the analytical 
form of I c i and I c 2 , 

= r<B 



e (A 0 — \E — eH|) 

\JE 2 - Aq^Aq — [E — eV) 2 



x Ci 



— cos 



X cos 



2 E-eV 
v F 

2 E-eV 

( 

v F 

2 hex 



d sin 



Si 



2 E-eV 



v F 

2 E-eV 



v F 



d - 



Ic 2 — 



d , 

u F J 

i r v/2 ~ A 



(B-7) 



dE 

0 (eV - 2A 0 ) 



^(E + eV/2? - A lj{E- eV/2 f - 
2 hex 



x cos — d 1 cos 

V v F 



V F 



-d , 



(B-8) 



where <tq = 167re 2 TgAL(0)AR(0) [mv/(27rd)] . 
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